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We study a Langevin equation derived from the Michaelis-Menten (MM) 
phenomenological scheme for catalysis accompanying a spontaneous repli- 
cation of molecules, -which may serve as a simple model of cell-mediated 
immune surveillance against cancer. We examine ho-w t-wo different and sta- 
tistically independent sources of noise - dichotomous multiplicative noise 
and additive Gaussian -white noise - influence the population's extinction 
time. This quantity is identified as the mean first passage time of the sys- 
tem across the zero population state. We observe the effects of resonant 
activation (RA) and noise-enhanced stability (NES) and -we report the ev- 
idence for competitive co-occurrence of both phenomena in a given regime 
of noise parameters. We discuss the statistics of first passage times in this 
regime and the role of different pseudo-potential profiles on the RA and 
NES phenomena. The RA/NES coexistence region brings an interesting 
interpretation for the gro-wth kinetics of cancer cells population, as the 
NES effect enhancing the stability of the tumoral state becomes strongly 
reduced by the RA phenomenon. 
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1. Introduction 



Since tlie times of Delbriick [Ij, a considerable researcli has been done 
on auto-catalytic reactions using stochastic models [2j-[4J. Spontaneous 
random molecular fluctuations play a significant role in the driving kinetic 
mechanism of many regulatory enzymatic reactions. Random fluctuations 
in genetic networks are inevitable as chemical reactions are probabilistic and 
many genes, RNAs and proteins are present in low numbers per cell. Gene 
expressions, for example, involves a series of single-molecule events subject 
to significant thermal fluctuations, which may result in a divergence of fate 
and lead to non-genetic population heterogeneity. Moreover, mean levels 
of heterogeneity for the populations are ensured by regulatory factors that 
tune the spontaneous fluctuations. Although many examples of potentially 
noise-exploiting cellular processes have been discussed and documented in 
literature, how cells function and process information when the underlying 
molecular events are random still remains an open question [3]-|14j. Other 
cellular processes influenced by noise include ion-channel gating, neural fir- 
ing, cytoskeleton dynamics and molecular motors [15j-[18j. 

Simulation schemes that model each individual reaction event can be 
computationally demanding [2] as the copy of the input and output molecules 
increases. Therefore, as an alternative way to account for molecular fluc- 
tuations in biochemical reactions is to assume that the time evolution of 
a species of interest is produced at an average rate /(x), degraded at rate 
g{x) and follows the stochastic differential equation: 



In the above equation the noise term is appended as the additive (white) 
noise source ^{t) weakly perturbing deterministic evolution of the system. 
Obviously, implicit for such a formulation is a continuous description of 
molecular species where the dynamics is cast in terms of infinitesimal changes 
of concentrations. Eq. ([1]) will be further postulated to study dynam- 
ics of a catalytic reaction . We consider the phenomenological Michaelis- 
Menten scheme for the catalysis accompanying a spontaneous replication of 
molecules 
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where a substrate X forms first a complex Z with molecules of the enzyme 
y, before the conversion of X to a product P is completed. By assuming 
that the production of X-type molecules inhibited by a hyperbolic activation 
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(Michaelis-Menten inhibition kinetics) is the slowest process under consider- 
ation and by considering a conserved mass of enzymes Y + Z = E = const, 
the resulting kinetics can be recast in the form of the 1-dim Langevin equa- 
tion 

with the pseudo- "free-energy" potential profile U{x) expressed as [I9] 

U{x) = -— + — + px-PHx + l), (4) 

where x is the normalized molecular density with respect to the maximum 
number of molecules, and with the following scaling relations 

x = —, = —, (3 = —, t = Xt. (5) 

The time-scale separation performed on a kinetic scheme of the system 
([2]) requires that a subset of species is asymptotically at steady state on 
the time scale of interest, and is known as a quasi-steady-state assumption 
[11|5] . It reduces the model complexity by effectively reducing the number 
of considered molecular species and reactions. 

In the case of Michaelis-Menten kinetics this approximation assumes also 
much larger concentration of substrate X than the enzyme concentration 
E and with the steady-state constraint dZ/dt = leads to a characteristic 
regulatory term —f3x/{l + x) in the dynamic ([3]). The resulting potential 
profile eq-dH) has at most three extrema representing deterministic station- 
ary states of the system (see Fig. [1]) 

XI = 0, (6) 

29 ' 



X. ^ (S) 



The essential feature captured by the model is, for a constant parameter 
9, the /3-dependent bistability. In the above form and by assuming time 
dependent, random variations of the parameter /3, the model has been used 
to describe an effect of cell-mediated immune surveillance against the cancer 
[20| . Most of tumoral cells bear antigens which are recognized as strange 
by the immune system. A response against these antigens may be mediated 



4 



acta'ania printed on February 4, 2008 




Fig. 1. The Michaelis-Menten potential with parameters: f3 — 3,9 — Q.l, A — 0.02. 
Labels: "o": extrema of C/(a;):x = 0, 4, 5; "x": extrema of C/+(a;): a; = 0, 4.28, 4.72; 
"*": extrema of C/-(x): x = 0, 3.83, 5.17. 



either by immune cells such as T- lymphocytes or other cells, not directly re- 
lated to the immune system (like macrophages or natural killer cells). The 
process of damage to tumor proceeds via infiltration of the latter by the 
specialized cells which subsequently develop a cytotoxic activity against the 
cancer cell-population. The series of cytotoxic reactions between the cyto- 
toxic cells and the tumor tissue may be considered to be well approximated 
by a saturating, enzymatic-like process whose time evolution equations are 
similar to the standard Michaelis-Menten kinetics. The variability of ki- 
netic parameters defining this process naturally affects the extinction of the 
tumor [20tl2T]. 

The process of population growth and decay can be described as a mo- 
tion of a fictitious particle in a potential U{x), whose shape varies in time 
due the noisy /? parameter. Minima of the potential correspond to stable 
states of the population, where, in case of a purely deterministic system, 
the cell concentration does not grow or decay. One of the minima is always 
at X = 0, which reflects the fact that the state in which cancer cells are 
absent is a stable one. The other minimum corresponds to a fixed-size cell 
population. Noise-induced transition between those stable states across the 
potential barrier may be interpreted as a "spontaneous" extinction of cancer 
due to environmental fluctuations, as well as due to random variations in 
the immune response efficiency. We will study a system with parameters 
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chosen in such a way that x = is a global minimum. We are interested 
in the lifetime of the metastable state. Mean escape time has been inten- 
sively studied in order to characterize the lifetime of metastable states of 
static and fluctuating potentials with different initial conditions [22]-[41]. 
The studies show that the mean escape time has different non-monotonic 
behaviors as a function of both the thermal noise intensity and the mean 
frequency of potential fluctuations. These behaviors are a signature of two 
noise-induced effects, namely the resonant activation (RA) [22j-[33j and the 
noise enhanced stability (NES) [S]-[3I]- The NES effect increases the aver- 
age lifetime of the metastable state, while the RA phenomenon minimizes 
this lifetime. 

The purpose of this work is to find the optimal range of parameters in 
which the positive role of resonant activation phenomenon, with respect to 
the cancer extinction, prevails over the negative role of NES, which enhances 
the stability of the tumoral state. The following section explains more 
extensively the concepts of resonant activation and noise-enhanced stability 
and their role in the analysis of cancer growth kinetics. In the next section 
we report the appearance of RA and NES phenomena in this model and we 
discuss the statistics of first passage times in this regime. In the same section 
the role of different potential profiles on the RA and NES effects, with piece- 
wise linear and MM potentials with reflecting boundary is analyzed. Finally 
we draw the conclusions. 



Our starting point is the following Langevin equation of an overdamped 
Brownian particle moving in a random fluctuating potential barrier 



Here ^{t) is a Gaussian process with zero mean and correlation function 
(^(*)?(*')) = (^^^{t - i')- The potential Vix,t) is the sum of two terms: 
the fixed potential U{x) and the randomly switching term G{x)r]{t), where 
r]{t) stands for a Markovian dichotomous noise switching between two levels 
{A"'",A~}, where = A and A~ = —A. The probability of switching 
(ibA — > =F^) per time unit, or the mean frequency of switching, is labelled 
u = l/(2r), and the autocorrelation function of the dichotomous noise yields 



2. The Model System 

2.1. Mean first passage time analysis 



dx 
'di. 




dx 




(9) 
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The two noise sources are statistically independent, i.e. (^(t)^/(s)) = 0. The 
potential V{x,t) therefore flips at random time between two configurations 

U^{x) = U{x) + G{x)A^. (10) 

The corresponding Fokker-Planck equations which describe the evolution of 
probability density of finding the state variable in a "position" x at time t 
are 



dtpix,A'^,t) = dx 



dx 2 



p(x,A±,t) 



- ^p(x,A±,t) + ^p(x,AT,t). (11) 
It It 

In the above equations time has dimension of [length]'^ /energy due to a 
friction constant that has been "absorbed" in a time variable. With the 
initial condition 

p{U±,Xs,t)\t=o = ^S{x-Xs), (12) 

from Eqs. (jlip we get the equations for the mean first passage times (MF- 
PTs): " "" ' ' 

_1 - ^ dU+{x)dT+{x) ^2 d^T+{x) 

T T dx dx dx"^ 

_^ ^ r+(x) _ T2(x)_ _ ^ JU-{x)dT-{x) ^ ^, d^T-{x) ^ ^^^^ 
T T dx dx dx"^ ' 

where T~^{x) and T~{x) denote MFPT for U^{x) and U~{x), respectively. 
The overall MFPT for the system is derived from the above set of equations 
T{x) = r+(x)+T~(x) by assuming a reflecting boundary condition at a; = a 
and an absorbing boundary at x = b: 

dTHx)_ _ 

, \x=a — U, 

dx 

T±(x)U=b = 0, (14) 



Analytical solutions of Eqs. (jl3p can be expressed in a compact form only for 
piece- wise linear or piece- wise constant potentials j22t I23j . More complex 
cases require either use of approximation schemes [261 [2S1 [29] , perturbative 
approach [27] , or direct numerical evaluation methods [251 Ell E2] ■ 



acta'ania printed on February 4, 2008 



7 



2.2. Resonant activation 



The signature of the resonant activation phenomenon is the occurrence 
of a minimum of the MFPT as a function of the mean frequency of the 
switching z/ of the multipUcative dichotomous noise. As a consequence the 
resonant activation effect minimizes the average lifetime of a population in 
the metastable state. Let us assume that the Brownian particle is behind 
the potential barrier, in a neighborhood of the metastable state. When the 
barrier fluctuations are very fast the particle "sees" the average potential 
between the higher and lower configurations. The MFPT will tend to a 
constant value corresponding to the average static potential 



If the barrier fluctuations are slower than the actual escape rate (r large), 
the particle will escape before any barrier flip occurs. Therefore, the mean 
first passage time also tends to a constant, which now will be an average of 
the escape times for the higher and lower configurations of the potential 



If the potential barrier is high enough or the additive noise is weak enough, 
the mean first passage time in a static potential can be approximated by 
the inverse of the Kramers escape rate, which increases exponentially with 
AU{x)/a'^, where AU{x) is the height of the barrier. Such an exponential 
dependence guarantees the existence of a minimum of MFPT (r). Its value 
at intermediate r is lower than both asymptotic mean first passage times, 
for very large or very small r [30]. Specifically, the MFPT at very low mean 
switching frequency (Eq. (|16|) ) will be higher than the MFPT at the high 
frequency limit (Eq. (|15|) ). In the middle frequency regime the effective 
escape rate over the fluctuating barrier is the average of the escape rates 



where K+ = 1/MFPT{U+) and = 1/MFPT{U~), with K_ <. K+. 
Because of the exponential dependence, the MFPT in this intermediate 
frequency regime will be smaller than both MFPTr-too and MFPTr^o- 

2.3. Noise-enhanced stability 

The signature of the noise-enhanced stability effect is the occurrence 
of a maximum of the MFPT as a function of the additive noise intensity, 



lim MFPT{U+, U~,t) = MFPT + 




(15) 





(17) 
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for a fixed mean frequency u of the multiplicative noise. As a consequence 
the NES effect, differently from the RA phenomenon, maximizes the aver- 
age lifetime of the population in the metastable state. The nonmonotonic 
behavior of the mean escape time as a function of the additive noise inten- 
sity a depends on the potential profile parameters, on the parameters of 
the multiplicative dichotomous noise and also on the initial position of the 
Brownian particle |35j-|38j. Noise enhances the stability of the metastable 
state with different peculiarities related to different dynamical regimes: the 
average lifetime can greatly increase when the noise intensity is very low 
with respect to the height of the barrier and the initial positions of the 
Brownian particles are in the "divergent" dynamical regime |4H . When 
the particle starts from initial positions within the potential well, at small 
values of a with respect the height of the potential barrier and for a fixed 
potential, it will follow the monotonic behavior of the Kramers formula, as 
a function of the noise intensity. The mean escape time will then tend to 
infinity for cr ^ 0. If the particle starts from outside the well, at small 
a it will at once run down from its starting region towards the absorbing 
barrier and its mean escape time will be approximately equal to the escape 
time of a deterministic particle. At high noise intensities, the mean escape 
time decreases monotonically, regardless of initial positions. At intermedi- 
ate noise intensities, a particle starting from the outer slope may sometimes 
be trapped into the well. Such events, although rare, can significantly in- 
crease the mean escape time because a trapped particle stays then in the 
well for a relatively long time. 



where x{t) is the concentration of the cancer cells. The profile of the cor- 
responding quasi-potential (Eq. [4]) presents a double well with one of the 
minima at x = 0. The region for x > can show either a monotonic be- 
havior or a local minimum, depending on the values of parameters 6 and 
f3. In the present investigation we used only parameters able to give a local 
minimum of the Michaelis-Menten potential for x > 0: = 0.1 and /3 = 3 
(see Fig. [1]). For x — > oo the potential shows a strong cubic repulsion. Our 
biological system is an open system whose random interactions with the en- 
vironment and the random fluctuations of the immune system are described 
by an additive noise term (J^{t) in Eq. (jlSp . and a dichotomous Markovian 
noise r/(t), of amplitude A and mean correlation time r in the (3 parameter. 



2.4- Population model 
The kinetics of our biological system is described by the equation 




(18) 
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As a consequence the stochastic Michaehs-Menten (MM) potential switches 
between two conformational states U^{x) 

U^i^) = -y + — + (/? ± A)(x - ln(x + 1)), (19) 
and the Langevin equations for the system are 

= 1(1 - to) -(li± A)-^ + «*). (20) 

X + 1 

The process of population growth and decay can be described as a motion 
of a fictitious particle in the switching potential between two configurations 
U^{x) and U~{x). For negligible additive noise and small concentration 
of cancerous cells, this model resembles a standard Verhulst equation with 
perturbing multiplicative dichotomous noise, which exhibits a complex sce- 
nario of noise-induced transitions, observable in a pattern of the stationary 
probability density [42j. Here, we will address kinetic properties of this 
model by studying the mean first passage time citebork between high and 
low population states in the system. We will study how the two different 
sources of noise as well as the position of the starting point Xin influence 
the mean first passage time. We put the absorbing boundary at x = and 
the reflecting one at x = oo. The event of passing through the absorbing 
boundary is equivalent to a total extinction of cancer. 



3. Results 

3.1. Michaelis-Menten and piece-wise linear potentials with absorbing 

boundary 

Here we analyze the role of different potential profiles with a reflecting 
barrier placed on the minimum of the fixed MM potential, that is in absence 
of dichotomous noise. We consider two cases: (i) MM potential, and (ii) 
piece-wise linear potential. Both cases are approximations of the effective 
MM potential (Fig. [TJ Particularly the piece-wise linear potential is a good 
theoretical model, which enable to obtain closed analytical expressions of 
the MFPT (see Refs. [22ii23j). In the following Fig. [2] we report the results 
of our simulations on both cases and compare the behaviors of MFPT for 
NFS and RA phenomena with those obtained using the exact MM potential 
profile (see Fig. [1]). The NFS effect, for a fixed value of the correlation 
time of the multiplicative noise and an unstable initial position, is shown 
in Fig. [2^. The reflecting barrier placed in the proximity of the metastable 



10 



acta'ania printed on February 4, 2008 



state pushes the particle away from it and therefore decreases the MFPT for 
all values of the noise intensity, producing an overall reduction of the NES 
effect. The Brownian particle experiences only two different slopes in the 
piece- wise linear potential, while more slopes are experienced, and therefore 
more time spent around the metastable state, in the MM potential [35j. This 
explain why the piece-wise linear potential is the worst approximation. For 
noise intensity values greater than the height of the potential barrier, the 
different shape of the potential profile becomes irrelevant and the decreasing 
of MFPT values is due to the reflecting boundary only. In Fig. [2)3 is reported 
the MFPT as a function of the correlation time of the multiplicative noise, at 
a fixed value of the noise intensity, for the same different potential profiles. 
The effect of the reflecting boundary is again to decrease all the MFPT 
values in the MM potential and therefore to reduce the absolute values of 
the RA phenomenon. The minimum of MFPT is of course more pronounced 
in the exact MM potential and the piece-wise linear potential is again the 
worst approximation. Both effects are robust enough to be observed in 
different potential profiles. 



3.2. Coexistence of noise- enhanced stability and resonant activation 

In this section we look for the coexistence region of RA and NES ef- 
fects. We calculate therefore the MFPT for the following initial positions 
Xin = 3, 3.65, 5. We performed a series of Monte Carlo simulations of the 
stochastic process (|2Up with an absorbing boundary at x = 0, natural re- 
flecting boundary at a; = +oo and the values of parameters: (3 = 3, 6 = 0.1, 
A = ±0.02. The statistics for each MFPT has been taken from = 10'^ (for 
Xin = 3.65) and A^ = 10^ (for Xin = 3 and 5) simulation runs. The results re- 
ported here are an extension of those presented in [33] . The results confirm 
the existence of resonant activation and noise-enhanced stability phenom- 
ena in the studied system. Moreover, we have shown that in a certain range 
of parameters both effects can occur together [381 IMl SO] . 

In Figs. [3] and m we present the combined view of RA and NES effects. 
The considerations of subsection 12. 21 are valid not only for trajectories start- 
ing from inside the potential well, but also for arbitrary initial positions in 
the potential, if only the particle has a chance to be trapped behind the 
potential barrier for some time. If the additive noise intensity is very large, 
a particle starting from the outer slope of the barrier does not "feel" the 
barrier at all. If, in turn, the additive noise is very weak, the particle slides 
down the slope in an almost deterministic way. But at intermediate values 
of a, the particle can (in some, very infrequent, realizations of the stochastic 
process) be trapped behind the barrier and, in case of such a rare event, 
its escape time changes non-monotonically as a function of r, in a way 
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Fig. 2. (a) MFPT as a function of the additive noise intensity cr^ for three potential 
profiles: [1] piece-wise linear with a reflecting boundary, with slopes ki — +0.32 and 
k2 = —0.025 in the average position; [2] Michaelis-Menten (MM) with a reflecting 
boundary; [3] exact MM potential. The unstable initial position is at xq = 3.4 and 
T = 0.1. (b) MFPT as a function of the correlation time of the multiplicative noise. 
The initial position is Xq — 0.1, cr — 0.1 and the number of realizations is = 10^. 
All the other potential parameters are the same as in Fig. [TJ 
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Fig. 3. MFPT as a function of the additive noise intensity <j and the correlation 
time of multipUcative noise. Initial positions: — 3 (left panel), Xin — 5 (right 
panel). Parameters: f3 = 3, 9 = 0.1, A = ±0.02. Number of simulation runs per 
each point: A'' = 10^. 




Fig. 4. MFPT as a function of the additive noise intensity a and the correlation 
time of multiplicative noise. Initial position: x-^ — 3.65. The parameters are the 
same as in Fig. [31 Number of simulation runs per each point: N = 10^. 



described in subsection 12. 2i This effect of coexistence of noise-enhanced 
stability and resonant activation effects can be observed in Fig. HI When 
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Fig. 5. Left panels: MFPT as a function of r for chosen values of a in the noise- 
enhanced stability regime. From top to bottom: a = 0.10, a = 0.08, a = 0.06. 
Right panels: Corresponding first passage time histograms. Colors from white to 
black denote a fraction of all trajectories. Black color denotes values from 0.1 to 
1. Initial position: = 3.65. The other parameters are the same as in Fig. [31 
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the trajectory starts on the outer slope of the potential, only small fraction 
of the realizations are trapped into the well, and thus the RA effect is very 
small because only few trajectories will contribute to it. In order to obtain 
a well-visible RA (a large enough number of trapped trajectories) we have 
to choose the initial position sufficiently close to the top of U~{x). We have 
chosen Xjn = 3.65. For this initial position we observe a co-occurence of 
NES and RA. Since both effects act here in an opposite way, there exists 
a regime of a and r parameters where noise-enhanced stability is strongly 
reduced by resonant activation. This region is well visible in Fig. |3]for val- 
ues of parameters: a ~ 0.06, r ~ 100. In Fig. [5] we present cross-sections 
of Fig. [3] at chosen values of a and histograms of the first passage times 
corresponding to each point on the graphs. Most trajectories (more than 
90% of all) are those which do not fall into the potential well. The trajecto- 
ries trapped are minor contibutions: their fraction (of order of 0.1, 0.01 and 
even 0.001) depends on the barrier switching rate and on the additive noise 
strength. At 0.01 < a < 1, the first passage times of particles, which had 
not been trapped into the well, are of order of 10. The first passage times 
of the trajectories trapped behind the potential barrier depend of course on 
r: In the region where the resonant activation is better visible, for small r 
the first passage times are of order of 10^, for large r they are of order of 
10^ and more, and for the intermediate values of r the first passage time 
drops down to the order of 10^. 



4. Conclusions 

We studied a Langevin equation derived from the phenomenological 
Michaelis-Menten scheme for catalysis accompanying a spontaneous replica- 
tion of molecules. It contains an additive noise term (Gaussian white noise) 
and a multiplicative noisy driving (dichotomous noise) in the term responsi- 
ble for inhibition of population growth. The models of that type are widely 
used in biophysical description of cell-mediated immune surveillance against 
cancer [19\ 1201 [2T\ H6| ITT] , analysis of cell survival including repair kinetics 
after exposure to ionizing radiation [45] or delineation of factors responsible 
for the efficient action of biological regulatory networks [31 El El E] • We 
examined how the two different sources of noise influence the population's 
extinction time, identified with the mean first passage time of the system 
over the zero population state. We observe appearance of two noise-induced 
resonant phenomena in the system: Mean first passage time displays here 
a non-monotonic behavior which depends on the characteristic parameters 
of the noises present in the system, namely the additive noise intensity a 
and the correlation time r of the multiplicative noise. The occurence of 
a minimum in MFPT (resonant activation) is related to the r parameter, 
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whereas the emergence of the maximum in MFPT (noise-enhanced stabihty) 
is connected with a. We report the evidence for co-occurrence of resonant 
activation and noise-enhanced stabihty in the same regime of noise param- 
eters. We observe that the strong stabihty enhancement of the metastable 
state comes from very few contributions from the trajectories trapped in 
it. Their first passage times are several orders of magnitude longer than 
the first passage times of the trajectories which approached the absorbing 
boundary without being trapped in the metastable state. We analyzed the 
role of different potential profiles, with a reflecting barrier, on RA and NES 
phenomena. We found a very interesting, from the point of view of growth 
kinetics of cancer cells, coexistence region of both effects. In this region the 
NES effect, which enhances the stability of the tumoral state, is strongly 
reduced by the RA effect, which enhances the cancer extinction. In other 
words, an asymptotic regression to the zero tumor size may be induced by 
controlling the noise affecting hyperbolic inhibition of a spontaneous prolif- 
eration of cells. 
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